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1.  INTRODUCTION 


The  end-of-life  operation  of  spacecraft  is  often  characterized  by  the  need  for  careful  electrical 
power  management  This  requires  accurate  knowledge  of  payload  electrical  requirements,  eclipse 
occurrence  and  duration,  battery  conditions,  and  solar-array  output  capability.  Typically,  the 
solar-array  output  has  degraded  in  a  steady,  well-understood  manner  determined  by  the  array 
components  and  the  radiation  properties  of  the  orbit.  This  study  was  undertaken  near  the  end  of 
life  of  a  vehicle  when  such  power  management  methods  were  in  use.  The  unusual  aspect  of  this 
system  was  the  observation  that  the  solar  array  produced  more  output  current  when  pointed 
somewhat  away  from  the  sun.  The  deviation  from  zero  sun  angle  that  produced  the  maximum 
output  was  empirically  found  to  depend  on  the  exact  bus  voltage  set  point  and  the  age  of  the 
array.  Since  the  array  was  continuing  to  age,  and  accurate  output  current  predictions  were 
needed,  a  study  was  commissioned  to  model  the  array  behavior. 
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2.  MODEL  DETAILS 


The  principal  aspects  of  the  solar-array  model  are  the  solar  cell  beginning-of-life  (BOL)  and  end- 
of-life  (EOL)  electrical  parameters,  and  the  steady-state  temperature  of  the  array.  The  important 
electrical  parameters  are  the  short-circuit  current,  Isc.  the  open  circuit  voltage,  V^,  and  the  cur¬ 
rent,  Imp,  and  voltage,  Vmp,  at  the  maximum  point.  The  temperature  coefficients  of  these 
parameters  and  the  effect  of  radiation  on  the  parameters  are  critical  as  well. 

2.1  ARRAY  TEMPERATURE 

The  temperature  of  a  solar  array  can  be  generally  described  by  the  function 


T  = 


(q-n)qcos8  +  q' 
EG 


(1) 


where  a  is  the  solar  cell  absorptivity,  tj  is  the  cell  efficiency,  q  is  the  solar  intensity,  0  is  the  sun 
angle,  q'  accounts  for  the  earth's  albedo,  e  is  the  total  array  emissivity,  and  a  is  the  Stefan- 
Boltzmann  constant  In  order  to  compute  the  temperature  directly,  each  of  these  terms  must  be 
known.  For  the  present  problem,  temperature  telemetry  data  were  available  from  sensors  located 
on  the  backside  of  the  array.  These  data  were  used  to  simplify  Eq.  (1)  to 


■tfKcosG’U 

1  +  K  J  ’ 


where  K  is  a  constant,  and  T0  is  the  temperature  at  zero  sun  angle.  Figure  1  shows  the  quality  of 
fit  to  the  flight  data  that  this  functional  form  provides.  This  agreement  was  judged  to  be  adequate 
for  the  problem.  Ground  test  results  and  a  detailed  thermal  model  had  shown  previously  that  a 
10°C  differential  exists  between  the  front  and  back  of  the  panel  under  a  wide  range  of  steady- 
state  and  transient  conditions.  Therefore,  the  function  in  Eq.  (2)  was  used  to  predict  the  backside 
temperature,  and  the  frontside  temperature  representing  the  solar  cells  was  assumed  to  be  10°C 
above  this. 
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Figure  1.  Amy  temperature  from  flight  data  (symbols)  and  model  (line). 


2J2  FORM  OF  ARRAY  IV  CURVE 

The  array  current-voltage  behavior  (IV  curve)  was  assumed  to  be  related  to  the  individual  solar 
cell  IV  curves  by  the  following. 


V array  —  Vceu  *  Nserjes  -  V<jjode  -  VjR 


(3) 


Iarray  —  Icell  *  (^parallel  *  Nfaii)*Isun*F*COS  0  , 


(4) 


where  the  N;  refer  to  the  number  of  cells  in  series,  in  parallel,  or  failed  open,  as  labeled;  Isun  is  the 
relative  solar  intensity;  and  F  is  a  factor  that  includes  miscellaneous  losses  such  as  cell  matching  at 
assembly,  adhesive  darkening,  etc. 

The  IV  curve  of  a  solar  cell  was  modeled  by  the  expression  (Ref.  1) 


I(V)  =  (Isc  -  B*(exp((Vcejj/C)  -  1)), 


(5) 


where 


V array  —  V^n  +  +  ViR 


(6) 
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C  =  (Vmp- Voei/lnd  -UpA^) 


(7) 


B  =  (Isc-ImP)*exp(-Vmp/C). 


(8) 


The  advantage  of  this  expression  is  that  it  only  requires  knowledge  of  the  three  prime  points  (I*, 
Voc,  Pmax)  of  the  cell  in  order  to  generate  the  complete  IV  curve.  It  has  been  found  to  accurately 
represent  the  behavior  of  Si  solar  cells  at  both  beginning  of  life  and  after  significant  radiation 
degradation,  and  over  an  appropriate  range  of  operating  temperatures. 

2 3  SOLAR  CELL  PRIME  POINTS 

In  order  to  reproduce  the  array  effects  observed  in  the  flight  data,  it  is  essential  to  have  accurate 
knowledge  of  the  solar-cell  prime  points  as  a  function  of  temperature  and  radiation  dose.  The 
temperature  dependence  is  customarily  expressed  in  terms  of  temperature  coefficients,  which  may 
vary  in  different  temperature  ranges.  For  a  constant  temperature  coefficient  k,  the  equation 
describing  the  variation  of  a  parameter  is  simply 


X(T)  =  X(28°  C)*k*(T  -  28°  C) , 


(9) 


where  X  is  the  prime  point  variable.  The  radiation  degradation  is  based  primarily  on  empirical 
data  and  typically  takes  the  form 


X/X0  =  l-A*ln(l-nJ>/<{)0) , 


(10) 


where  A  and  «t>0  are  parameters,  X0  is  the  beginning-of-life  value  and  <j>  is  the  actual  equivalent 
1  MeV  electron  fluence. 

This  rather  simple  view  of  these  effects  is  complicated  by  the  second-order  effect  of  temperature 
coefficient  variation  with  radiation.  However,  it  is  found  that  at  the  high  radiation  doses  this  array 
has  experienced  (on  the  order  of  6E 15  equivalent  1  MeV  electrons/cm2),  the  temperature  coeffi¬ 
cients  have  essentially  reached  asymptotic  with  respect  to  radiation  effects. 

The  values  of  the  parameters  in  Eqs.  (9)  and  (10)  are  obtained  from  qualification  tests  performed 
on  the  solar  cells  prior  to  their  use.  For  the  K6  3/4  silicon  solar  cells  used  on  the  subject  array, 
there  were  two  qualification  test  documents  available.  These  are  referred  to  as  the  Preflight  Test 
and  the  Hughes  Miniqual  Test  (Ref.  2).  Table  1  compares  these  parameters. 

Given  the  differences  in  these  values  and  some  questionable  temperature  behavior  of  the  Preflight 
Test  data,  the  Miniqual  data  were  selected  for  use.  A  10%  uncertainty  was  assigned  to  all  of  the 
temperature  coefficients. 
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Table  1.  Comparison  of  K6  3/4  Solar  Cell  Parameters  from  Two  Qualification  Tests 


Parameter 

Preflight 

Mini-qual 

Temp.  Coeff. 

Isc  (A/cm2) 

6.50E-05 

5.90E-05 

Voc(V) 

-2.13E-03 

-2.25E-03 

^mp  (A/cm2) 

2.70  E-05 

4.20E-05 

Vmp(V) 

-2.41  E-03 

-2.22E-03 

Radiation 

Uc 

A 

0.0593 

0.0543 

<f>0 

5.25E+13 

3.41  E+13 

< 

o 

o 

A 

0.0312 

0.0272 

*0 

1.20E+13 

1.15E+13 

Imp 

A 

0.107 

0.0582 

4»o 

2.54E+14 

4.18E+13 

Q. 

E 

> 

A 

0.0306 

0.0258 

4*o 

9.63E+12 

7.79E+12 
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3.  FLIGHT  DATA 


The  flight  data  shown  in  Table  2  were  provided  for  validation  of  the  model.  The  output  current 
of  the  +Y  solar  array  wing  was  measured  as  a  function  of  solar  incidence  angle  and  bus  voltage. 
The  data  were  taken  in  rapid  succession  to  minimize  changes  in  sun  angle  and  albedo.  The  range 
specified  for  each  angle/voltage  combination  is  the  95%  confidence  interval  based  on  multiple 
measurements.  The  uncertainty  in  the  bus  voltage  values  is  150  mV. 


Table  2.  Current  (amps)  from  +Y  array  wing. 


Bus  voltage 

30.13  V 

30.53  V 

31.06  V 

31.53  V 

Angle 

15 

19.78-19.86 

19.30-19.40 

18.51-18.59 

17.76-18.06 

18 

19.79-19.89 

19.37-19.45 

18.57-18.69 

17.98-18.11 

23 

19.73-19.83 

19.26-19.43 

18.63-18.95 

17.96-18.14 

27 

19.56-19.76 

19.27-19.39 

18.73-18.81 

18.19-18.41 

35 

18.93-19.03 

18.48-18.64 

18.07-18.21 

4.  ANALYSIS 


At  the  time  the  flight  data  were  taken,  the  equivalent  1  MeV  electron  fluence  received  by  the  solar 
array  was  6.24E15  cm  2  for  Isc  and  1.35E16  cnr2  for  Voc,  lmp,  and  Vmp.  Figure  2  shows  the  IV 
curve  for  the  array  at  BOL  (dashed  curves)  and  at  present  (solid  curves)  for  two  temperatures. 
These  temperatures,  43°C  and  S4°C,  are  the  frontside  array  temperatures  at  solar  incidence  angles 
of  0°  and  30°,  respectively. 

The  parameters  in  Table  1  are  not  listed  explicitly  in  the  qualification  test  documents.  Rather, 
they  have  been  obtained  by  fitting  curves  to  graphical  representations  of  the  data  in  the  docu¬ 
ments.  The  temperature  coefficients  in  the  table  are  valid  only  in  the  range  0  <  T  <  60°C  and 
for  radiation  doses  (equivalent  1  MeV  electron  fluence)  >  1E15  cm*2.  Both  of  these  constraints 
are  satisfied  in  the  present  problem.  The  radiation  degradation  parameters  A  and  $0  are  obtained 
by  fitting  Eq.  (9)  to  the  data. 

The  observed  angle  dependence  of  the  array  output  can  be  understood  qualitatively  as  follows. 
The  power  system  is  shunt  regulated  and  designed  to  operate  at  approximately  30  to  33  V.  This 
voltage  range  operates  the  solar  cells  to  the  left  of  Pmax.  Figure  2  shows  no  significant  tempera¬ 
ture  dependence  of  output  current  in  this  region  at  BOL.  However,  at  the  high  radiation  doses 
experienced,  these  operating  voltages  lie  near  or  to  the  right  of  Pmax.  The  output  current  of  the 
degraded  array  is  then  very  sensitive  to  temperature  because  of  the  proximity  of  the  operating 
voltages  to  Voc.  Figure  2  illustrates  the  large  variation  in  current  that  can  be  caused  by  a  tempera- 


Figure  2.  Broken  lines  show  current-voltage  behavior  for  array  at  two  tempera¬ 
tures  at  beginning  of  life.  Solid  lines  are  for  present  state  of  degradation. 
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ture  change  under  these  operating  conditions.  Therefore,  under  some  conditions,  it  might  be 
expected  that  when  the  array  is  pointed  off  from  the  sun  vector,  the  gain  in  current  resulting  from 
the  lower  array  temperature  might  overcome  the  cosine  factor  loss  in  current. 

Figure  3  shows  the  flight  data  in  Table  2  (symbols)  and  the  predictions  of  the  model  detailed 
above  (curves).  Two  symbols  are  used  for  each  data  point,  representing  the  upper  ai._  rower  95% 
confidence  limits  shown  in  the  table.  All  of  the  data  were  measured  in  the  same  week  in  1992 
with  the  exception  of  the  32.6-V  data  set  which  was  taken  in  1991.  The  agreement  between  the 
model  and  the  data  is  quite  good  in  two  important  respects.  First  the  sun  angle  that  produces  the 
maximum  current  output  is  correct  for  all  three  bus  voltages.  Second,  the  magnitude  of  the 
maximum  is  well  predicted.  These  are  the  two  most  important  parameters  used  in  determining 
the  maximum  array  output  capability  and  how  to  obtain  it.  Given  this  information,  a  low-risk 
power  management  plan  can  be  implemented  for  use  in  eclipse  season. 


Figure  3.  Model  predictions  (curves)  and  flight  data  (symbols)  for  off-pointed 
array  performance  as  a  function  of  sun  angle.  Double  symbols  are  used  to  indi 
cate  95%  confidence  limits  of  flight  data. 
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5.  SUMMARY 


It  has  been  observed  that  a  solar  array  that  has  suffered  significant  radiation  degradation  pro¬ 
duces  maximum  current  output  at  non-zero  sun  incidence  angles.  An  analysis  has  been  pre¬ 
sented  that  shows  that  this  effect  can  be  accurately  modeled  using  radiation  degradation  and  tem¬ 
perature  coefficients.  The  net  gain  in  current  can  be  as  great  as  25%  under  realistic  operating 
conditions.  Furthermore,  knowledge  of  the  array  performance  characteristics  enables  a  low-risk 
power  management  plan  to  be  formulated  for  use  during  eclipse  season. 
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